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i , 1 INTRODUCTION 



ABSTRACT 

We employ numerical simulations of galaxy mergers to explore the effect of galaxy mass 
ratio on merger-driven starbursts. Our numerical simulations include radiative cooling 
of gas, star formation, and stellar feedback to follow the interaction and merger of four 
disk galaxies. The galaxy models span a factor of 23 in total mass and are designed to 
be representative of typical galaxies in the local Universe. We find that the merger- 
driven star formation is a strong function of merger mass ratio, with very little, if 
any, induced star formation for large mass ratio mergers. We define a burst efficiency 
that is useful to characterize the merger-driven star formation and test that it is 
insensitive to uncertainties in the feedback parameterization. In accord with previous 
work we find that the burst efficiency depends on the structure of the primary galaxy. 
In particular, the presence of a massive stellar bulge stabilizes the disk and suppresses 
merger-driven star formation for large mass ratio mergers. Direct, co-planar merging 
orbits produce the largest tidal disturbance and yield that most intense burst of star 
formation. Contrary to naive expectations, a more compact distribution of gas or 
an increased gas fraction both decrease the burst efficiency. Owing to the efficient 
feedback model and the newer version of SPH employed here, the burst efficiencies of 
the mergers presented here are smaller than in previous studies. 
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While mergers between galaxies of approximately equal 
mass, so called major mergers, garner a significant amount 
of observational and theoretical interest, it is likely that 
galaxies accreting much smaller objects are also important 
drivers of galaxy evolution. This possibility is a natural con- 
sequence of the hierarchical growth of structure predicted 
by the Lambda cold dark matter (ACDM; a.k.a., "Double 
Dark") cosmology (see, e.g.. ISomerville et al.ll20oil ). In this 
scenario every dark matter halo, and presumably the galaxy 
it hosts, is assembled though a continuous and varied merger 
history comprised of m any small accretion events and a few 
major" ones (see, e.g.. lLacev fc Colelll993l : IWechsler et al.l 



whose shapes are determined by random motions, i.e., they 
are pressure-supported, and which possess projected mass 
distributions that scale approx i mately as 



-1/4 



(e.g. , 

Barnes fc Hcrnquist Il99ll. IT9 92; Hcr nquistJ 1 1991 Il993bl: 



Naab fc Burkerd 120031 ; ICox et all l2006l l?h. Observations of 



mai 
2002; 



I)- 

The transformation of galactic systems by major merg- 
ers has been well studied. Numerical simulations demon- 
strate that major mergers between spiral galaxies pro- 
duce remnants that generally resemble elliptical galaxies, 
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fine structure, kinematic subsystems, the bimodal distribu- 
tions of global clusters in elliptical galaxies, and the light 
profiles of merger remnants lend support to the assertion 
that elliptical galaxies are the byproduct of major galaxy 
mergers - in what is co mmonly termed the "merger hy- 
pothesis" (|Toomrel 119771 ) . For further information regard- 
ing the observational and theoretical evidence supporting 
the merger hypothesis, or the general dynamics of galaxy 
interactions, we refe r the i nterested reader t o revi ews by 
iBarnes fc Hernquistl i|l992l ), ISchweizerl l| 19981 ), and IStruckl 
( 20061 ). 



While major mergers are likely to be the most dramatic 
events in a galaxy's assembly history, it is likely that the 
more prevalent minor mergers give rise to the wide vari- 
ety of g alaxy morpho l ogy that defines many classification 
systems |Hubblelll926l : Ide Vaucouleursll 19771 : Ivan den Berghl 
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Il96d ; ISandaed Il975l ). The morphological transformation 
owing to the effects of minor merger has been studied 
via numerous numerical simulations, which demonstrate 
that galactic disks become warped a nd heated as a re- 
sult of the accretion of a satellite (IQuinn fc Goodman 



198d : lToth fc Ostrikerll99l:lQuinn et al.lll993l:IWalker et al 



1996 



Hua ng fe Carlberd 1 19971 ; IVelazquez fc White! [l999: 



Font et al.l 2001), that the satellite itself is tidally stripped 
to produce stel l ar st r eams throughou t the galactic halo 



to produce stel l ar st r eams t hroughout the galactic halo 
llJohnston et all Il999l; iHelmi fc White! l200ll ; llVlaver et al.1 
|2002| ; iBovlan-Kolchin fc Mai I2007T ). and that the rem- 
nant galaxy will be a systematically earlier Hubble type 
ilAguerri et al.ll200ll;lNaab fc Burkertl|2003l; iBournaud et alj 
120041 , 120051 ; lEliche-Moral et al.1 120061 ; iMaller et all l200rj ). 
Furthermore, there is an increasing amount of evi- 
dence that these theoretical expectations ar e consistent 
with observations in the loca l univ e rse (e.g.. Ilbata et al.l 
200ll; iDalcanton fc Bernstein! |2002| ; iNavarro et all ' 



Bullock fc JohnstonlbOOsI ) 



While both major and minor mergers are expected 
to produce some degree of morphological evolution, it is 
still unclear whether all mergers will lead to a merger- 
driven starburst and therefore drive spectrophotometric evo- 
lution. There is mounting evidence that near equal mass 
(major) mergers trigger the most vigorous star-forming 
galaxies in the local Universe, the ultraluminous infrared 
galaxies (ULIRGs) with bolometric lumino s ities greater 
10 12 L Q (see, e.g., ISanders fc Mirabell Il996l ; iBorne et al.1 
2000). More generally, large statistical samples consistently 
indicate a clear anticorrelation between projected galaxy 
aair separation and star formation ind i cators in the optical 



pair s eparation and star lormation indicato rs m the o ptical 
llBarton et al.ll2q00l; [Nikolic et al.ll2004l ; |Lambas et al.ll2003l: 
Sol Alonso et al.1 120061 ; iBarton et al ' 



2007) and the infrared |Geller et al.l l200rj ; ISmith et al 



2007: Ellison et al 



20071 ). Moreover, numerical simulations naturally explain 



these observations as centrally concentrated starbursts trig- 
gere d by tidal forces that attend the galaxy interaction (see , 
e.g., jBarnes fc Hernquistl Il99ll; iMihos fc Hernquistl Il994b . 



19961; ISpringelll2000l; iKapferer et al.l 120051 : ICox et al.1 12006 ; 
di Matteo et al.ll2007l ). 

In contrast, relatively little is known about the re- 
sulting starbursts during interactions between galaxies of 
unequal mass. Observationally, statistical studies are ham- 
pered by magnitude limits and confusion with background 
sources, but preliminarily indicate that enhanced star for- 
mation (as measured by Ha) is strongly dependent upon 
the relative mass (or magnitude) of the interacting galax- 
ies. Initial results using the CfA2 Redshift Survey indicated 
no correlation be tween galaxy pair separation and equiva- 
lent width of Ha (| Woods et al.ll2006l ), suggesting that minor 
mergers do not lead to merger-driven star formation. More 
recent studies using the Sloan Digital Sk y Survey (SDSS, 
lEllison et ai]|2007l ; IWoods fc Gelleil 120071 ) find that minor 
mergers indeed lead to enhanced star formation, however 
the data also suggest that the enhancement is primarily 
observed in the satellite galaxy while the larger primary 
galaxy has star forming properties simil ar to non-interacting 
field galaxies (|Woods fc Gelled 120071 ). More locally, de- 
tailed observations of individual galaxies in the nearby Uni- 
verse suggest that minor mergers may be responsible for 
some degree of star for mation enhan cement, e.g. , NGC 
278 (|Knapen et al.l I2004T ), NGC 3310 (|Smith et all 1 19961 ). 



NGC 4064 and NGC 4424 (|Cortes et al.ll2006l). NGC 7742 
jMazzuca et al.1 120061), and the M 81 group jWalter et al.1 
120021 ; iForster Schreiber et al.ll2003l ). 

While numerical simulations have been used extensively 
in the study of star formation in major mergers, there exist 
only a handful of previous simulations that have specifically 
quantified the star formation history of unequal mass galaxy 
mergers. These works have established that tidal perturba- 
tions can induce bar formation, inflows of gas, and subse- 
quent bursts of star formation when the primary disks are 
susceptible to i nstabilities, such as when they do not contain 
a stellar bulge j|Hernquistlll989l ; IMihos fc Hernquistl [l994al ; 
iHernquist fc Mihosl I1995F T However, because these studies 
performed only a limited number of simulations, and did not 
include some important physical processes, our understand- 
ing of the relationship between mergers of unequal mass 
galaxies and the resulting star formation is still highly in- 
complete. 

In this paper we fill in this gap by describing a large set 
of numerical simulations of unequal mass disk-galaxy merg- 
ers. In parti cular, we seek to extend and improve the pioneer- 
ing work of IMihos fc Hernquistl (Il994al ) and provide a more 
complete picture of merger-driven star formation, including 
its dependence upon the merging orbit and the structure 
of the primary galaxy. We will also take advantage of re- 
cent improvements in the treatment of supernovae feedback 
(|Springelll200"ol ; ISpringel fc Hernquisdl2003l ; ICox et al.ll2006l ) 
and smoothed-particle hydr odynamics (SPH) methodology 
ijSpringel fc Hernquistl I2002D which allow us to follow the 
evolution of galaxies with larger gas fraction than was pos- 
sible previously. 

The results of our simulations are relevant for a num- 
ber of related studies. Many models of galaxy formation of- 
ten employ merger-driven star formation as a necessary in- 
gredient to produce bright galaxies throughou t cosmic time 
jGuiderdoni et al.lll998l ; ISomerville et al.ll200ll ; iBaugh et al.l 
120051 ). Our results will serve as a useful input to such models. 
In addition, our models will be a valuable tool for interpret- 
ing observations. By comp arison to the observ e d star forma- 
tion i n pairs of galaxies (IBarton et al.l |2000| ; IWoods et al.l 
120061 ; IWoods fc Gelled [20071 ; iLin et al.M2007l ). we can con- 
strain our models and infer the contribution of mergers 
to the universal star formation rate. The numerical mod- 
els have also been used to generate "simulated observa- 
tions" of the merging galaxies through radiative transfe r 
modeling (| Jonssorfeooj ; I Jonsson et alj|2006l ; |jonssonll2006l). 
to st udy dust attenuation in spiral galaxies ( Rocha et al.l 
120071 ). to understan d the size, shape, and scaling relations 
of merger remnants ijPekel fc Cox] I2OO6I ; iNovak et all 120061 ; 
ICovington et al.1l2007l ). and to calibrate non-parametric in- 
dicators of galaxy morphology (Lotz et al., in preparation), 
which can then be used to quantify the merger rate over 
cosmic time. 

We organize this paper as follows: Our numerical tech- 
niques are summarized in SjS] an d the galaxy models used 
for our galaxy interactions are described in £j3j A general 
description of the mergers is provided in i|4] including a dis- 
cussion about the possible methods to quantify the merger- 
driven star formation. Some additional models are explored 
in ij5] Finally, we discuss our results and conclude in £|6] 
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2 NUMERICAL SIMULATIONS 

All of the numerical simulations per formed in this wor k 
use the N-Body/SPH code GADGET (ISpringel et al.ll200ll) . 



Coxl l|2004 ) and 



In general , pur methodology is similar to 
ICox et all (|2006h . and thus, in this section we will only re- 
view selected aspects of our techniques and methodology 
that are relevant to this work. We use the first version of 
GADGET, however the smoothed particle hydrodynamics 
(SPH) modules are upgraded (with the aid of V. Springel) 
to use the "conservative entro py" version that is detailed 
in ISpringel fc Hernq uist (2002]). The radiative cooling rate 
A ne t(p, u) is comput ed for a primordial plasma as described 
in lKatz et al.l fl996). Stars are formed at a rate determined 
by the local SPH density normali zed to match observed star 
formation rates l|Kennicutj| 19981 ) . Furthermore, we employ 
a threshold density pth, below which stars do not form. Fi- 
nally, star formation occurs within individual SPH particles 
and each particle can stochastically spawn two stellar par- 
ticles. 

One of the largest uncertainties associated with current 
numerical simulations including star formation is the im- 
plementation of supernova feedback. Because of the limited 
resolution achievable in state-of-the-art numerical simula- 
tions, most work performed to date adopts a "sub-grid" 
approach where the physical processes associated with feed- 
back are included in a simple, yet flexible manner with free 
parameters that can be tuned to match observations. Unfor- 
tunately, there are often a range of acceptable free parame- 
ters which produce distinct st ar formation histories or rem- 
nant morphologies (see, e.g.. [Thacker fc Couchmanl |2000l ; 
iKav et alj|200a ISpringel et alj|2005bl ; ICox et alj|2006f ). For 
many of the simulations performed in this work, we have run 
two different feedback models, either the nOmed or n2med 
models as introduced in lCox et al.l l|2006T ). Both models are 
termed "medium" because they dissipate feedback energy 
on a 8 Myr timescale. The n2med model treats star-forming 
gas with a stiff equation of state while nOmed assumes this 
gas is isothermal with an effective temperature ~ 10 5 K. We 
will compare these two feedback models in £14.3.21 in order 
to determine how these assumptions influence the resulting 
star-formation history. 

All simulations presented throughout this work use a 
gravitational softening length of 0.1 kpc for all baryonic 
particles, and 0.4 kpc for dark matter particles. The SPH 
smoothing length is required to be greater than half the 
gravitational softening length, or > 50 pc. We have per- 
formed a few tests with both smaller and larger soften- 
ing lengths and found little or no differences in the star- 
formation history, sugge sting that we have achieved numer- 
ical convergence (see lCox et al.ll2006l . for similar tests). 



3 GALAXY MODELS 

Because our primary motivation is to investigate collisions 
between unequal mass galaxies, we wish to construct galax- 
ies with a range of different masses. Furthermore, because 
the properties of galactic disks depend on mass, we also 
need to design galaxy models that are physically plausible. 
To begin, we select the largest galaxy model to have a total 
stellar mass of M* ~ 5.0xl0 10 Mq, slightly above the tran- 
sitional mass of 3.0xl0 10 M Q found to divide large galaxies 



with old stellar populat ions from smaller star forming ones 
|Kauffmann et akll2003l ). Three smaller galaxy models are 
included that have M* equal to 1.5 xlO 10 Mq, 0.5 xlO 10 
Mq and 0.1 xlO 10 Mq, making the total stellar mass span 
a range of 50. We label the largest of these models G3, and 
the rest G2, Gl, and, finally, the smallest GO. 

Once the stellar masses are fixed, we next determine 
the stellar size using the late-typ e galaxy size-mass relation- 
ship found by IShen et al.l (120031 . eq. 18), from the analysis 
of 140,000 galaxies in the Sloan Digital Sky Survey. This re- 
lation fixes the half-light radius R50 for each model, which 
we assume is equivalent to the stellar half-mass radius. 

Because each model is assumed to be a late-type galaxy, 
the y consist o f a stel lar and gaseous exponential disk, a stel- 
lar iHernquistl (|l990ft bulge, all embedded in a dark matter 
halo. A general description of the methods employed to con- 
struct model dis k galaxies, i nclud ing the assumed profiles, 
can be found in ICox et al.l ([20061 ) ■ however, those models 
are not identical to the ones used here. The model parame- 
ters used for this work, which are listed in Table [T] are se- 
lected such th at each mo del galaxy is statistically average, 
as opposed to lCox et alj (|2006l ) who specifically modeled a 
gas-rich Sbc galaxy. 

To ensure that the models are statistically average, the 
mass of the stellar disk and bulge are constrained by the 
half-mass radius R50 and the observed b ulge-to-disk ratio 
of local late- type galaxies (|de Jondll996l ). To quantify the 
latter constraint, we assume a K-band mass-to-light ratio of 
0.7 an d 0.5 for t he bu lge and disk, respectively, and fit a line 
to the lde Jond (| 19961 ) data, finding 



log(M bu i gc ) = 1.6 log(M disk ) - 1.03, 



(1) 



where the masses are in units of 10 10 Mq. Using this formula, 
along with the fixed total stellar mass for each model, both 
the stellar bulge and stellar disk masses are determined. To 
determine the size of these components (the exponential disk 
scale radius R<j and the bulge scale radius Kb) an iterative 
approach is used where the first step is to guess an initial 
disk scale radius R^. The bulge scale radius R(, is then fixed 
by the e mpirical correl ation between disk and bulge sizes 
found bv lde J oriel (1 19961 . see his Eq. 5). This initial guess for 
the disk and bulge mass distributions is then compared to 
the desired half-mass radius R50 and corrected until the two 
are within 1% of each other. 

In addition to the stellar component, each galaxy 
model contains gas, the fuel for star formation. Moti- 
vated by observations of l ocal late-type galaxies (e.g., 
Bro eils fc van Woerdedfl994l) . gas is distributed in an ex- 
tended exponential disk with a scale radius proportional to 
that of the stellar disk, r gas = ar sta ,. In our default model we 
assume a — 3. To determine the mass of gas, we use the ob- 
served trend that sys tems with lower total/stellar mass have 
a higher gas f raction dBell et al.ll2003l ; lGeha et alll2006l '). Us- 
ing data from lBell et al.1 ( 20031 . kindly provided in electronic 
form by E. Bell), we parameterize this relationship as 



log(M g 



0.78 log(M*) + 1.74, 



(2) 



where both masses are in units of 10 10 Mq . Because this re- 
lation is derived from the mean gas content as a function of 
total baryonic mass, the model disks will contain only mod- 
erate quantities of gas, and are therefore unlikely to produce 
large bursts of star formation. In i]5.4l we will perform a small 
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Table 1. Properties of galaxy models used in this study. The compound galaxy has a total mass M v j r , and is composed of dark matter 
and a stellar disk, stellar bulge, and gaseous disk. The total mass of each baryonic component and their mass fractions are provided. The 
dark halo has an NFW profile with concentration c and spin parameter A. The stellar disk has a total mass which is a fraction 

m^ of the total mass and a fraction j d of the total angular momentum. The disk has an exponential profile of scale radius R d and scale 
height zq. Each galaxy model contains a spherical bulge of mass Mb u i gc and mass fraction mj,, with an exponential profile of scale length 
R 6 and zero angular momentum. The gas mass is M gas , a fraction m s of the total mass, and j s is the angular momentum fraction. The 
gas disk, like its stellar companion, is assumed to have an exponential surface density with a scale radius a times R^. We define the gas 
fraction / as the fraction of the total disk mass which is gaseous. Finally, each system is represented by N particles, which represent the 
dark matter, stellar disk, bulge and gaseous galaxy. 





G3 


G2 


Gl 


GO 


M* (10 10 M ) 


5.0 


1.5 


0.5 


0.1 


R50 (kpc) 


3.86 


2.88 


2.35 


1.83 


Total Mass, M vir (10 10 M ) 


116.0 


51.0 


20.0 


5.0 


Concentration, c=R v j r /r' s 


6 


9 


12 


14 


Spin Parameter, A 


u.uu 


U.LM 


u.uo 


U.UJ 


Disk Mass, M disk (10 10 M ) 


4.11 


1.35 


0.47 


0.098 


Disk Mass Fraction, mjj 


0.035 


0.026 


0.024 


0.019 


Disk Scale Length, R^ (kpc) 


2.85 


1.91 


1.48 


1.12 


Disk Scale Height, zo (kpc) 


0.40 


0.38 


0.30 


0.22 


I 1 1 Q w ^ T\T T\ H T f '^ 1 /~\T1 1 1 


0.015 


0.010 


0.010 


0.010 


Bulge Mass, M bu i g e (10 9 M e ) 


8.9 


1.5 


0.3 


0.02 


Bulge-to-disk ratio, B/D 


0.22 


0.11 


0.06 


0.02 


Bulge Mass Fraction, m b 


0.008 


0.003 


0.002 


<0.001 


Bulge Radial Scale Length, R;, (kpc) 


0.37 


0.26 


0.20 


0.15 


Gas Mass, M gas (10 10 M©) 


1.22 


0.48 


0.20 


0.06 


Gas Mass Fraction, m s 


0.011 


0.009 


0.010 


0.012 


Gas Fraction, f 


0.196 


0.242 


0.286 


0.375 


Gas Scale Multiplier, a 


3.0 


3.0 


3.0 


3.0 


Gas Spin Fraction, ] g 


0.012 


0.010 


0.013 


0.019 


N 


240,000 


150,000 


95,000 


51,000 


N dm 


120,000 


80,000 


50,000 


30,000 


N gas 


50,000 


30,000 


20,000 


10,000 


Ndisk 


50,000 


30,000 


20,000 


10,000 


Nbul g c 


10,000 


10,000 


5,000 


1,000 



number of additional mergers using galaxies with higher gas 
fractions. 

The final component of the galactic system is the mas- 
sive d ark-matter halo. Eac h model contains a dark halo 
with a lNavarro et alj (1 19961 , NFW) profile whose properties 
are selected so that the r otation curve satisfies the bary- 
onic Tully-Fisher relation (|Bell fc de Jond l200ll ; [Bell et all 
120031 ). Because most of the models require very little non- 
baryonic mass in order to satisfy the Tully-Fisher relation, 
we do not include adiabatic contraction. Still, the resultant 
halo concentrations c are b elow the mean found in cosmo- 
logical N-body simulations (|Bullock et al.ll200ll). a tension 
that has been emphasized previously (see. e.g.. lAlam et al.l 
2002). We note that (motivated by observations) the galaxy 
models have systematically higher baryon fractions in higher 
mass galaxies, resulting in a total mass ratio 23 between G3 
and GO as opposed to a ratio of 50 in stellar mass. 

The c o nstruc tion of these galaxy models is based upon 
iHernquistJ <|l993al ) and its more recent incarnation to in- 
clude a NFW halo l|Springel fc Whitelll999l ; ISpringelll2000t ). 
We note that because the spatial distribution of baryons is 
fixed by the galaxy model, the halo spin does not determine 
the size of the gaseous and stellar disk, as in many popular 



models for disk formation (see, e.g.. iMo et al.|[l99cf ). How- 
ever, we feel that our procedure is adequate since it pro- 
duces stable disks, provided that the velocity dispersion of 
the disks is chosen such that the Toomre stability parameter 
Q is greater than 1 at all radii and the feedback model can 
provide adequate pressure support to the gas. 

3.1 Evolution of Isolated Systems 

The stability of our model galaxies is demonstrated in Fig- 
ure[T] which shows the projected gas and stellar surface den- 
sities after being evolved in isolation for ~ 1 Gyr (greater 
than 5 orbital periods in all models). While the mass dis- 
tributions of each component are stable during the entire 
period that they are simulated, the noisy potential resulting 
fro m the finite partic le number seeds transient instabilities 
(see IHernquistJ Il993al ) . These instabilities are most promi- 
nent in the G3 and G2 galaxies, which have larger disk sur- 
face densities than the less massive galaxy models. 

The stability of the galaxy models is also reflected by a 
constant star-formation rate (SFR), regardless of the feed- 
back model employed, as shown in Figure [2] The average 
SFR is 0.95, 0.25, 0.06 and 0.001 M yr _1 for models G3, 
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Figure 1. Projected mass density for each galaxy model, GO (left column), Gl, G2, and G3 (right column) when simulated in isolation 
for ~ 1 Gyr, which is greater than 5 orbital periods at the half-mass radius in all models. Each row shows the projected stellar or 
gaseous mass as specified in the lower-left of each panel. The top two rows show the galaxies face— on, while the bottom two rows show 
the edge-on view. Each panel measures 80 kpc on a side. A color bar in the upper— left panel indicates the stellar surface density scale 
that is used for all panels, and the color bar in the one below indicates the gaseous surface density scale that is used in all panels. 



G2, Gl, and GO, respectively, during the 1 Gyr period shown 
in Figure [2] The average SFR scales as M ga s , as i s expected 
from our Schmidt-type star-formation law (see ICox et all 
l200fj ). The SFR of the smaller galaxies, Gl and GO, dis- 
play increased fluctuations owing to the more significant im- 
pact of feedback, decreasing resolution, and the fixed star- 
formation threshold density p t h- The SFRs shown in Fig- 
ure [2] are constant for ~ 2 Gyr. After this, galaxies G3 and 
G2 begin to consume a significant fraction of their high- 
density (> pth) gas and their SFR eventually drops to ap- 
proximately one-tenth its initial value at 6 Gyr. Another 
consequence of the adopted star-formation law is that the 
star- formation timescale scales with gas density, viz., tsf ~ 
Mg as /SFR ~ p~ 1//2 , in other words, more massive galax- 
ies, with higher gas densities, have shorter gas consump- 
tion timescales. As a result, the more massive galaxy models 
consume a larger fraction of their initial gas mass than the 
smaller ones. During the entire 6 Gyr evolution, model G3 
converts 26% of its original gas to stars, G2 converts 22%, 
Gl converts 11%, and GO converts 3%. 



4 GALAXY MERGERS 

The previous section defined four galaxy models G3-G0 and 
established the stability of these models by simulating them 
in isolation. In this section we discuss the outcome of binary 
mergers between these models on a single merging orbit. 
Performing a single merger between all combinations of the 
four models yields four equal-mass major mergers and six 
unequal-mass mergers. In the following section ([}5]), we dis- 
cuss a number of additional merger simulations to determine 
systematic dependencies. 



4.1 Fiducial Encounters 

For reference, we list the merger mass ratios in Tabled Mass 
ratios are given for the total, stellar and baryonic masses. 
In what follows we will label a simulated mergers by its 
two merging constituents. For example, G3G1 is a merger 
between our largest galaxy G3 and the second-smallest Gl. 
The total mass ratio is 5.8 to 1, i.e., model Gl is 5.8 times 
smaller than G3. When referring to the interacting galaxies, 
the less massive galaxy will be termed the satellite and the 
more massive galaxy the primary. We note that some galaxy 
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Figure 2. Star formation for our isolated galaxy models G3, G2, 
Gl and GO. Each model has been run with two separate feedback 
models and both star-formation histories are shown, nOmed with 
a dotted line and n2med with a solid line. 



Table 2. Total, stellar, and baryonic mass ratios between our 
four galaxy models G3, G2, Gl, and GO. 



Primary 


Satellite 


Total 


Stellar 


Baryonic 


G3 


G3 


1:1 


1:1 


1:1 


G3 


G2 


2.3:1 


3.3:1 


3.1:1 


G3 


Gl 


5.8:1 


10.0:1 


8.9:1 


G3 


GO 


22.7:1 


50.0:1 


38.9:1 


G2 


G2 


1:1 


1:1 


1:1 


G2 


Gl 


2.6:1 


3.0:1 


2.8:1 


G2 


GO 


10.0:1 


15.0:1 


12.4:1 


Gl 


Gl 


1:1 


1:1 


1:1 


Gl 


GO 


3.9:1 


5.0:1 


4.4:1 


GO 


GO 


1:1 


1:1 


1:1 



models (G2, Gl, and GO) can be a primary in one interaction 
and a satellite in another. 

All galaxy interactions are initialized with a nearly un- 
bound elliptical orbit with eccentricity e =0.95. The initial 
separation is fixed to be slightly less than the virial radius 
of the primary and is 250, 100, 80, and 50 kpc for G3, G2, 
Gl, and GO models, respectively. The pericentric distance is 
set to be 13.6, 3.8, 2.96, and 2.24 when the primary is G3, 
G2, Gl, and GO, respectively. The resulting interactions are 
fast and nearly radial, consistent with orbit s found for dark 
matter halos in cosmological simulations (IVitvitska et al 



20021 : iKhochfar fc Burker t]|2006l : lBenson1l200"5l :l Zentner et al 
20051 ). The only exception to the above orbits is the merger 
G3G0, where the total mass ratio is 22.7 to 1, the largest 
mass ratio we simulate. In this case, the fiducial interaction 
was still not fully merged after the system was evolved for 
~ 12 Gyr. In order to better compare this interaction with 



the other, more rapid, mergers, we multiplied the initial ve- 
locity of the satellite by 0.2 so that the merger occurred 
during the 6 Gyrs we followed each interaction. 

Finally, the fiducial series of orbits are all prograde, i.e., 
the angular momentum of the orbit is (nearly) aligned with 
the spin of the primary. A slight tilt of 30° is introduced such 
that the interaction is not entirely co-planar. In section [5TT1 
we resimulate the G3G2 and G3G1 mergers with a variety 
of orbits to determine the dependence of merger-driven star 
formation on orbital angular momentum (R pC ri) and merger 
alignment. 

4.2 Merger Evolution 

The evolution of a typical merger event is shown in Fig- 
ures [3] and [4] for the G3G2 merger, and Figures [5] and [6] for 
the G3G1 merger. In both cases, the projected stellar surface 
density is first shown, followed by the projected gaseous sur- 
face density. The images shown in these figures are typical of 
the remaining mergers that are not shown, and are also con- 
sistent with what has been discussed in studies of equal mass 
merg e rs (e.g.. | Barnes fc Hernauistlll99ll : lMihos fc Hernauistl 
1 19961 : ISpringell l2000l: ICox . et all l2006h. and unequal mas s 



mergers (e.g., iHernquistJ 1 19891 : iHernquist fc Mihosl ll9Q5h 
For completeness, the following provides a brief outline of 
the merger process, and the resulting star formation. 

In both of the mergers shown the satellite galaxy (G2 
and Gl) first makes a fast, direct approach toward the pri- 
mary galaxy (G3). This close passage tidally distorts the 
disks and generates symmetric tails in both. Owing to the 
energetic orbit, the satellite separates again for several or- 
bital periods (>1 Gyr) before returning for its second en- 
counter with the disk of the primary. A generic feature of all 
the interactions simulated here is the efficient loss of angu- 
lar momentum by the satellite. In particular, even for large 
mass ratio encounters the orbits become almost entirely ra- 
dial after the first or second close passage. 

The loss of orbital angular momentum and eventual ra- 
dial nature of the final merger limits the coupling between 
orbital and spin angular momentum, and therefore the tidal 
response. Most of the induced star formation during the fi- 
nal coalescence is a result of the collisional nature of the gas 
and the rapid fluctuations in the potential as the satellite 
oscillates through the center and the system finally relaxes. 
The collisional nature of the gas also produces several other 
features distinct from the stellar disk. During the first close 
passage, the satellite creates a snow-plow effect, as it clears 
out gas from the primary's extended gas disk. This effect 
continues during each close passage until eventually the ma- 
jority of the gas in the satellite has been stripped off, leaving 
a predominantly stellar satellite. However, by this point the 
orbit has decreased such that it is difficult to disentangle the 
two galaxies. 

While it is beyond the scope of the current paper to 
provide a detailed analysis of the structure of the merger 
remnants, there are several qualitative features present in 
Figures [3H5] that motivate future inquiry. First, it is evident 
that the stellar disks have been thickened and dynamically 
heated during the accretion event, con sistent with a num- 
ber of previous numerical studies (e.g.jQuinn 



Goodman 



1986 |: lQuinn et al.lll993l:lToth fc Ostrikerll992l:IWalker et al. 
19961 : IVelazquez fc White! Il999l : iBenson et all 12004 
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Figure 3. Projected stellar mass density during the G3G2 merger simulation as viewed in the orbital plane. Each panel measures 200 kpc 
on a side and the time, in Gyr, is displayed in the upper left of each panel. The orbit of the satellite galaxy G2 is denoted by a dotted 
line until it has completely merged with the primary galaxy. The bottom-right panel shows a side view of the final merger remnant, and 
clearly indicates the initial 30° inclination of the progenitor disk. The top— left panel indicates the color scale used in all panels. 
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Figure 4. Similar to Figure [3] but for the gaseous component. 
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Figure 5. Similar to Figure [3] but the G3G1 merger is now shown. 
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Figure 6. Similar to Figure [6] but for the gaseous component. 
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iBournaud etafl 12004 I2005T ). The side views demonstrate 
that the G3G2 interaction has had a more disruptive effect 
on the primary's stellar and gaseous disk than the G3Gf 
interaction, an unsurprising result given the smaller mass 
ratio and hence larger tidal disturbance. Furthermore, the 
gaseous disks are also affected. It is particularly interesting 
that the G3Gf remnant contains a relatively large, thin 
gaseous disk, indicating that during an unequal mass 
interaction the gaseous disk can dissipate its energy and 
maintain its structure while the stellar disk cannot. The 
second qualitatively interesting feature is the efficient 
deposition of the satellite galaxy in the center of the merger 
remnant. Such a process may be a p rominent mechanism 
to build galactic bulg es (see also lAguerri et al.l l200ll ; 
lEliche-Moral et aJjl200a ). 

4.3 Star Formation 

The star formation history for all of our fiducial galaxy merg- 
ers is displayed in Figure [7] There are four plots in this 
figure, one for each of the four primary galaxy models. The 
upper-left panel shows the star-formation rate (SFR) for all 
interactions that involve the most massive galaxy G3. The 
remaining panels show all the remaining mergers, grouped 
by the largest participant galaxy. Thus each panels shows 
the SFR for one major merger and one or more minor merg- 
ers, except for the GO panel which only has the GOGO major 
merger. 

The SFRs shown in Figure [7] are clearly enhanced com- 
pared to the quiescent SFRs shown in Figure For ex- 
ample, the maximum SFR during the G3G3 interaction 
(~ 26 Mq yr" 1 ) is 13 times larger than the summed SFR 
of two isolated G3 disks (~ I Mq yr" 1 ), the maximum SFR 
during the G3G2 interaction (~ 27 Mq yr" 1 ) is 21 times 
larger than an isolated G3 (~ 1 Mq yr" 1 ) plus an isolated 
G2 (~ 0.25 M yr" 1 ), and the maximum SFR during the 
GOGO interaction (~ 3.8 Mq yr -1 ) is 3800 times larger than 
two isolated GO disks (~ I0" 3 M yr" 1 ), however 34.3.21 
demonstrates that these factors depend upon the feedback 
model. 

In 34.41 we will quantify the relationship between the 
mass ratio of the interacting galaxies and the enhanced star 
formation in more detail, but we note that this trend fits nat- 
urally within a merger-driven scenario for star formation. In 
particular, during an interaction gas is stripped of its angu- 
lar momentum by bar-like structures during the early stages 
of the merger or the abundant collisions and a growing po- 
tential well during the messy coalescence. In both cases the 
resulting central concentration of gas fuels a burst of star 
formation. Since the tidal forces associated with the merger 
generate these effects, their magnitude is expected to scale 
with the size of the perturbation. While this general pic- 
ture has been studied in great detail for collisions between 
equal mass ga l axies (e.g 



Mihos fc Hernouistlll994al . 1 19961 ; 



ISpringell[2OO0l ; [c ox et al.ll2006r ). we demonstrate here that 
the star formation clearly depends on the mass ratio of the 
interacting galaxies. 

We caution that the global SFR, as is shown in Fig- 
ure [7] will be dominated by the largest galaxy when the 
interacting galaxy mass ratio is large, as in the G3G1 and 
G3G0 mergers. In this case, the small enhancement in the 
SFR is only an indication that the primary, G3, is not tidally 



P^h 0.6 



0.4 



0.2 
0.0 



G3G1, 1:5.8 merger 

— Total (both G3 and Gl) 
contribution from G3 




contribution from Gl 



/ 







2 3 4 
Time (Gyr) 



Figure 8. Star— formation rate during G3G1 interaction including 
the contributions from participant galaxies G3 and Gl separately. 



disturbed enough to induce radial inflows of gas and addi- 
tional star formation. The less massive satellite, however, 
is expected to experience significant tidal forces and may 
have a large enhancement of star formation that is hidden 
within the global rate. Such is the case for G3G1, as shown 
in Figure [8] where the global SFR is shown as well as the 
contributions from both the primary and the satellite. While 
G3, the larger galaxy, constitutes the majority of the overall 
star formation, its rate is largely unchanged from the initial 
stages and eventually decays as a result of gas depletion. 
There is a small enhancement of star formation above the 
quiescent G3 from T w 3 — 5 Gyr coincident with the final 
coalescence. In contrast, the satellite galaxy Gl experiences 
several bursts of star formation that are > 3 times its qui- 
escent rate, consistent with merger-driven star formation. 

The exact relationship between merger mass ratio and 
enhanced star formation is complicated by the fact that the 
star formation history also depends upon the nature of the 
interaction, the structure of the participating galaxies, and 
the "sub-grid" model for star formation and feedback. The 
first two dependencies are inherently physical assumptions 
and will be investigated further in Sj5] The last dependency, 
however, is determined by numerical free parameters that 
are not well constrained. In 34.3.21 this issue will be con- 
sidered further, with the specific goal of determining which 
features of the star-formation history provide a measure of 
the merger-driven star formation that are insensitive to as- 
sumptions about the feedback model. 



4-3.1 Location of Starburst 

One result of merger-driven star formation is a concentra- 
tion of newly formed stars in the galactic center. Figure [9] 
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Figure 7. Star formation for our fiducial minor merger series. The total mass ratio is listed next to each run label. The line thickness 
increases with decreasing mass ratio, e.g., the thickest line in each panel is the 1:1 major merger, while the thinnest line is the large mass 
ratio minor merger. Note that the vertical axis changes from panel to panel, based upon the maximum star formation rate for the major 
merger. All runs shown here used the nOmed feedback model. 



demonstrates that this is a generic outcome of our interac- 
tions by showing the surface mass density of "new stars," 
i.e., stars that form during the interaction, for the remnants 
of the G3Gx set of interactions. Numerous numerical simu- 
lations have previously found similar profiles (though often 
plotted versus r 1 ^) for r emnants of major mergers (e.g ., 
iMihos fc Hernquisti Il99rj ; ISpringell l200d : ICox et al l 120061) . 
The central concentrations of new stars are often disjoint 
from the outer, old-star profiles, in apparent conflict with 
the r 1 / 4 profile that most ellipticals are assumed to have. 
While this was once considered a problem for the merger 
hypothesis, recent high-resolution observations now indi- 
cate that such features exist in nearly all merger remnants 



|Rothberg fc Joseph! 120041) and man y low-luminosity ellip- 
tical galaxies jKormendv et al.ll2007f ). 

Because the gaseous inflows that produce these cen- 
tral concentrations of new stars are produced by the grav- 
itational forces arising from the interaction, we expect the 
amount of new stellar mass in the galactic center to correlate 
with the merger mass ratio. Indeed, Figure [9] demonstrates 
that the new star profile within ~ 4 kpc becomes monoton- 
ically steeper with decreasing merger mass ratio. The least 
disruptive mergers, G3G1 and G3G0, have profiles that are 
essentially exponential beyond ~ 1 kpc, while the smaller 
mass ratio mergers G3G3 and G3G2 have profiles that are 
better fit by r 1//4 -type distributions. 

Another way to characterize the distribution of star for- 
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Figure 9. Mass surface density of "new stars," i.e., stars that are 
forming during the interaction, for mergers G3G3, G3G2, G3G1, 
and G3G0. 



mation as a result of the interaction is by simply quantifying 
the fraction of gas that is driven to within the central sev- 
eral hundred pc of the primary. This quantity was an espe- 
cially relevant measure of the potential starburst event when 
numerical simu l ations did not include st ar formation (e.g., 
Hcrnquist 1989; lHerncmist fc~ Mihos 199^). In these works it 
was always found that > 35% of the gas in the primary was 
driven to within the central several hun d red p c. The only 
exception was when iHernquist fc Mihosl (|l995l ) adopted a 
model in which the interstellar medium was isothermal with 
a temperature > 3 x 10 K. We note that we always find 
that < 35% of the gas in the primary is driven to within 
500 pc, even during the G3G3 major merger. For the large 
mass ratio mergers this percentage is closer to 10%, which 
is similar to the G3 galaxy evolved in isolation. This trend 
toward smaller inflows of gas is another manifestation of the 
differences between our modeling and what was performed 
previously (see more extensive discussion in i )6.1[) . 

While it is beyond the scope of the present work, analy- 
sis of the remnant profiles uncovers two interesting avenues 
for future study. First, it is intriguing that the new-star 
profiles of the remnants of very large mass-ratio mergers are 
well fitted by a bulge plus exponential disk profile, suggest- 
ing that such minor mergers may be an efficient mechanism 
for growing galactic bulges. Second, at large (~ 10 kpc) radii 
the profiles of the large mass ratio remnants differ from the 
G3 evolved in isolation owing to induced star formation and 
angular momentum transfer during the early stages of the 
interaction (see also lYounger et alll2007t ). 



Dependence on Star Formation and Feedback 
Parameters 



4.3.2 



As me ntioned pr e vious ly (SjU and ij4.3p . and shown explic- 
itly bv lCox et ail (|2006l) using numerical simulations of ma- 
jor mergers, the star-formation history during an interac- 
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Figure 11. Comparison of the gas consumption during equivalent 
simulations performed with two different parameterizations for 
feedback, nOmed and n2med. The label for the merger simulations 
are grouped by the primary. The isolated galaxies are identical to 
the simulations presented in Figures [T] and [2] 



tion depends upon the adopted feedback model. This is in 
contrast to the star-formation history of an isolated sp iral, 
which, as shown in Figure [2] (and in ICox et al1l2006l ). is 
independent of the feedback model. As a result, it is unclear 
what measure of the star-formation history will provide a 
robust characterization of the merger-driven star formation. 
To investigate this further, we resimulated our fiducial set of 
mergers using the n2med feedback model in addi tion to the 
nOme d that was presented in £14-31 As shown in ICox et all 
(2006), the "medium" feedback models are slightly favored 
because they can maintain a stable gas disk in an isolated 
gas-rich Sbc galaxy model and also produce merger-driven 
starbursts. The nOmed model treats star-forming gas with 
an isothermal (~ 10 5 K) equation of state while the nlmed 
has a stiff equation of state that restricts the quantity of gas 
at very high densities. 

Figure [TU] shows the SFR and gas consumption during 
the G3G2, G3G1, and G3G0 unequal mass mergers for both 
the nOmed and n2med feedback models. For all interactions 
the peak of star formation begins sooner, has a much larger 
maximum, and a shorter duration in the nOmed feedback 
model as opposed to nlmed. However, regardless of the dif- 
ferences in the star-formation history, the gas consumption 
is very similar, as shown for the full series of merger simu- 
lations in Figure 1111 This result motivates us to focus our 
analysis upon quantities that use the gas consumption rather 
than details of the time-dependent star-formation history, 
since the gas consumption is invariant to assumptions about 
the feedback model. We also note that the non-equivalence 
in gas consumption between the two feedback models that is 
present in Figure [11] can be used to estimate our errors ow- 
ing to adopting a single feedback model. From the standard 
deviation in discrepancies we estimate an error of ~ 0.04 in 
the gas consumption owing to the choice of feedback model. 
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Figure 10. Comparison between two parameterization of star formation and feedback in identical unequal mass mergers. The top panels 
show the star-formation rate, and the bottom panels show the gas consumption. In general, details of the star-formation history depend 
on the feedback model while the gas consumption remains relatively invariant. The burst efficiency (see i|4.4,ll l of each run [nlraed on 
top, and nOmed on the bottom) is provided in the lower panel. 



However, we also caution that 3 out of the 4 major mergers 
have slightly (< 0.1) higher gas consumption when adopting 
the nOmed feedback model, and the 2 minor mergers with 
the largest mass ratios have higher gas consumption when 
adopting the n2med feedback model. Therefore, it is possible 
that subtle, yet systematic, trends exist between alternate 
feedback models. 

As a final comment we note that Figure \W\ demon- 
strates that the feedback model can drastically affect the 
time-dependent star-formation rate, and therefore the lu- 
minosity of the system, even though the total gas consump- 
tion is similar. The varying luminosity evolution admits the 
possibility to better constrain the feedback model through 
a comparison t o the observed distr ibution of star-formation 
rates (see, e.g., iNoeske & aU l20O7h. 



4-3.3 The Star-Formation Timescale 

The previous section demonstrates that each feedback 
model yields a unique star formation history, and noted 
that one difference between the histories is the star- 
formation timescale. Because differences in the star- 
formation timescale directly affect the predicted color, line- 
strength, and nuclear activity of the interacting pair, fu- 
ture work may have the ability to discriminate between 
feedb ack models based upon comparisons wit h observations 



ack m odels base d upon co mparison s wit l 
iBarton et all l2000t rNikolic et al.1 12004 



Woods et al 



Ellison et al 



20061; IGeller et alj|2006l; I Woods fc Gellerll2007l ; 
20071 ; barton et alj|2007r ). 

As a preliminary step toward this goal, we quantify 
the starburst timescale for the mergers simulated here from 
three different fits to the starburst events. The first fit is 
motivated by the intrinsic shape of the starburst events as 
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Figure 12. Star-formation rate during the G3G3 major merger 
with the nOmed feedback parameters (thick solid line). Overplot- 
ted are three characterizations of the merger timescale, one, a 
single Gaussian fit to the merger-induced peak of star forma- 
tion (thin solid line), two, a decaying exponential fit to the star- 
formation rate after the peak of star formation and, three, the 
full width at half maximum (vertical dashed lines). 



exemplified in Figures [7] and 1101 namely they appear Gaus- 
sian in shape and their timescale can be characterized by 
the width a. The second fit is motivated by many models 
of galaxy evolution that describe star formation as a decay- 
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Figure 13. The star formation timescale (tsf) as a function of 
merger mass ratio for our fiducial series of mergers when sim- 
ulated with both the nOmed and n2med feedback models. The 
timescale is the average of the three fits (Gaussian, exponential, 
and FWHM) scaled to the exponential decay timescale r and 
the errors are estimated by the dispersion between the three fit 
parameters. The solid lines are fits to the simulated data using 
Equation [3] 



ing exponential with a timescale r, and are therefore often 

Somcrvillc fc Primackl 



described as "tau" -models (see, e 



1 19991 ; lHarker et all |2006| ; iNoeske fc all \200H) . The last 



method to describe the star formation timescale calculates 
the full width at half maximum (FWHM) and therefore 
avoids any biases associated with adopting a specified dis- 
tribution for of the starburst. Another possibility, which we 
reserve for future work, is to compile the full luminosity dis- 
tri bution for the e ntire merger as has been done previously 
bv lJonssonl l|2004 ) and lHopkins et all l|200fj ). 

As an example of a starburst event, and its correspond- 
ing fits, we show the merger G3G3 in Figure [12] Printed di- 
rectly on the figure are estimates for the starburst timescale 
that result from each fit. The three separate timescale mea- 
sures can be placed on an equal footing by noting that 
the FWHM of a Gaussian distribution is 2.35a, and fitting 
an decaying exponential to one side of a Gaussian yields 
r ~ l.lcr. For this particular case, the starburst has a dis- 
tinct peak and extended tails compared to the Gaussian fit 
(i.e., a positive Kurtosis). Therefore, the timescale derived 
from the FWHM is much smaller than that derived from 
the Gaussian fit. The starburst is very well-fitted by the 
tau-model (so much so that the fit is difficult to see in Fig- 
ure I12[) which yields a slightly longer estimate of the star 
formation timescale. 

In general, we note that nearly all models that employ 
the nOmed feedback have peaky star formation, and thus 
small FWHMs compared to the Gaussian or exponential 
estimates. However, because the star formation events are 
typically irregular, we consider the best estimate of the star 
formation timescale to be an average between r, l.lcr, and 



0.47FWHM, i.e., the estimates from the different fits are 
all scaled to be in terms of the exponential decay timescale 
t. This also provides a simple estimate of the error in our 
fitting procedure from the dispersion among the three de- 
terminations of the timescale. 

In Figure[Jj|]we plot the estimates for the star formation 
timescale and the associated errors for the fiducial set of runs 
as a function of merger mass ratio for both feedback models. 
To characterize the mass-ratio dependence, we have fit a 
simple power-law to the mergers of each feedback model, 



tsf 



to 



Ms. 



M 

nrin 



(3) 



Confirming our visual inspection of the SFRs shown in Fig- 
ure [10] the n2med model has a much longer star-formation 
timescale than the nOmed model. In additional, both models 
display a trend for the star-formation timescale to increase 
with increasing mass ratio, however the uncertainties are 
large enough that this trend is not definitive. In particular, 
the goodness of fit is only slightly reduced if the fit is fixed 
to have no mass-ratio dependence at all (a = 0). 



4.4 Parameterizing Star Formation Enhancements 
During Galaxy Interactions 

The star- formation histories shown in Figure [7] demonstrate 
clear signs of enhanced star formation during the merger. 
Furthermore, i|4.3.2l argues that this enhanced star forma- 
tion is more robustly characterized by the amount of gas 
consumed, rather than time-dependent quantities such as 
the maximum star formation rate, owing to uncertainties in 
the adopted feedback model. However, we must also con- 
sider that the two galaxies which participate in the interac- 
tion would have converted some of their available gas into 
stars even in isolation, and therefore the gas consumption 
itself does not provide a complete characterization of the 
merger-driven star formation. 

A more explicit comparison between the merger-driven 
star formation and that of the isolated disks is presented 
in Figure 1141 which shows several measures of star forma- 
tion for the G3G3, G3G2, G3G1, and G3G0 mergers, in 
each column, from left to right, respectively. The top row 
presents the star-formation rate (SFR) for the merger as 
well as the summed SFR of the primary and satellite galax- 
ies when evolved in isolation. For clarity, the difference be- 
tween the two SFRs is shaded. The second row shows the 
ratio of the merger SFR to the combined SFR of the iso- 
lated galaxies, clearly delineating when the rate is enhanced 
owing to the merger. The third row plots both SFRs divided 
by the stellar mass, which is also the inverse of the stellar- 
mass doubling time, and is sometimes termed the "specific 
SFR." The fourth row plots the SFRs divided by the gas 
mass, which is the inverse of the gas consumption timescale. 
Finally, the bottom row presents the gas consumption of 
both the merger and the combined isolated galaxies with 
the difference shaded. 

There are several features present in Figure [14] that de- 
serve more discussion. First, reinforcing the notion first pre- 
sented in Figure [7] there is a clear correlation between inter- 
acting galaxy mass ratio and merger-driven star formation. 
Significant episodes of star formation only occur for mergers 
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Figure 14. The rows from top to bottom show the star-formation rate (SFR), SFR normalized by the combined SFR of both isolated 
galaxies, the SFR divided by the total stellar mass, the SFR divided by the total gas mass, and the gas consumption for mergers 
G3G3, G3G2, G3G1, and G3G0, from left to right. Panels show quantities during the interaction with a solid (red) line and for the 
combined isolated galaxies with a dotted (green) line. The difference between the two lines is shaded to emphasize the merger-driven 
star formation. All simulations shown employ the nOmed feedback model. The burst efficiency (see EI4.4.H is printed in the bottom row 
for each interaction. 



between galaxies of nearly equal mass. For mergers between 
galaxies with a large mass ratio, such as G3G0, the cumula- 
tive star-formation history is nearly indistinguishable from 
the primary G3 evolved in isolation. 

By comparing the star-formation history of G3G3 to 
that of G3G2, we uncover another interesting feature of Fig- 
ure 1141 Namely, these two interactions have similar peak 
levels of star formation, even though the total amount of 
gas consumed is much less for the higher mass ratio G3G2 
interaction. This feature of the star-formation history em- 
phasizes that caution needs to be exercised when attempting 
to quantify the merger-driven star formation. 

As a final comment, we note that all models except 
G3G0 show elevated star formation at the end of the inter- 



action (T ~ 6 Gyr) compared to the primary G3 evolved in 
isolation, as indicated by the normalized SFR shown in the 
second row of Figure 1141 Therefore, most measures of the 
the merger-driven star formation are subject to uncertain- 
ties that depend upon the duration over which the simula- 
tions follow the merger (if the simulation was run for long 
enough, all galaxies would consume essentially all of their 
initial allotment of gas). In practice, these errors are quite 
small, which we estimate to be ~ 0.03 from the differential 
gas consumption (bottom row of Figure I14p if the simula- 
tion is followed for 1 Gyr prior to, or later than, the current 
6 Gyr we adopt as a standard. 
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where e\-\ is the burst efhciency for equal mass mergers, 



M sat and M D 



are the mass of the satellite and primary, 
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Figure 15. Burst efficiency e (defined in Eq. U) for all fiducial 
merger simulations. Interactions in which G3 was the largest pro- 
genitor, i.e., G3Gx, are shown with a (red) asterisk, while G2Gx 
mergers are (green) open triangles, GIGx are (blue) X's, and 
G0G0 is an open (black) square. The solid line is the best— fit re- 
lation using Eq.[5] and the dotted line is the best-fit relation using 
Eq. [6] The vertical line below the key shows our estimated error 
in burst efficiency of ±0.07 owing to uncertainties in the feedback 
model and the duration in which the simulation is followed. 



4-4-1 Burst Efficiency 

In order to quantify the merger-driven star formation, we 
wish to define a simple, useful, and robust parameterization 
of the star formation owing solely to the interaction. To 
this end, we are motivated by the discussion of ^4.3.21 to 
introduce the "burst efficiency" e, as 



(fraction of gas consumed during interaction) 
(fraction of gas consumed by constituent 
galaxies evolved in isolation, during the 
same time period). 



(4) 



The burst efficiency is the fractional amount of stars that 
are formed (or equivalently, the gas that is consumed) dur- 
ing the merger that would not have formed in isolation. 
Figure [T5] shows the burst efficiency for our fiducial merger 
series and reaffirms what was implicit in Figures [71 and [141 
namely, merger-driven star formation is only significant dur- 
ing mergers where the participating galaxies are comparable 
in mass. 

The burst efficiency appears to be a smoothly decreas- 
ing function of the merger mass ratio and its parameter- 
ization may be useful for future studies of galaxy forma- 
tion. In fact, such a parameterization has alread y been in- 
troduced by ISomerville. Primack fc Faberl (|200ll . see their 
Fig. 1) based upon a much smalle r set o f data (the simula- 
tions of lMihos fc HernquistlH994al.ll996l), and used in their 
semi-analytic models. Following ISomerville et al.1 (|200ll ). we 
adopt the following form for the burst efficiency: 



Ms 



respectively, and 7 fixes the mass ratio dependence. Perform- 
ing a least-squares fit to our entire set of fiducial mergers we 
find ei : i=0.55 and 7=0.69. This fit is shown as a solid line 
in Figure [Pol 

Owing to the relatively small star-formation enhance- 
ment in large mass ratio mergers (M sat /M pr i mary < 0.2), 
and the systematic overestimate in the regime provided by 
the best-fit Eq. [5] we also consider an alternate form for 
the burst efficiency that is fixed to zero when the mass ratio 
M sat /Mprimary is below eo and is 



M s 



Mprin 



eo 



(6) 



e = ei:i 



v M p r in 



(5) 



when the mass ratio M sa t /M pr imary is greater than eo- Per- 
forming a least squares fit to our fiducial mergers yields 
e 1:1 =0.56, 7=0.50, and eo=0.09, although the reduced x 2 
is nearly unchanged if eo is manually set to anything less 
than 0.11 (including 0, in which case the fit is identical to 
Eq. [5}. In other words, our fiducial set of mergers are con- 
sistent with there being no enhancement of star formation 
below a mass ratio of 9:1, although the data doesn't neces- 
sarily require this. The best fit to Eq.[6]is shown in Figure [Tol 
as a dotted line. 

A number of other fitting formula are also possible for 
the burst efficiency, such as single and broken power-laws, 
and various polynomials, however with the small number of 
data points and the associated uncertainties in calculating 
the burst efficiency, no one formula was statistically better 
than any other formula. Hence, we adopt Eq. [5] to describe 
the burst efficiency. 

For comparison, the ISomerville et al.l (|200ll ) fit to the 
Mihos & Hernquist data using Eq. [5] yielded era = 0.75 
and 7 = 0.18, when the progenitor disks did not include 
a stellar bulge, and era = 0.75 and 7 = 1.16 when they 
did. The consistency of e\-\ indicates that the presence of 
a bulge does not affect the burst efficiency of major merg- 
ers. We also note that a value of 0.75 for era is about 50% 
larger than our preferred value, indicating that the models 
employed here are less efficient at turning gas into stars than 
those of Mihos & Hernquist. In fact, this differ ence is similar 
to that found previously by ICox et all l|2006l ) who showed 
that the feedback model and the newer version of SPH em- 
ployed produce these ch anges. Our value of 7 is between the 
ISomerville et al.1 i|200ll ) value, a trend that is likely driven 
by the fact that our primary galaxy G3 has a bulge-to-disk 
ratio that is lower than in their model. We will attempt to 
explore these dependencies in more detail in the following 
section. 



5 ADDITIONAL MERGER SIMULATIONS 

The previous section detailed the outcome of 10 mergers, 
including 6 unequal mass mergers, between our four galaxy 
models. In this section we perform a number of additional 
merger simulations in order to explore several interesting, 
and particularly relevant regions of parameter space; the 
merging orbit and the structure of the progenitor disk, 
specifically the extent of the gaseous disk, the bulge-to-disk 
ratio, and the gas fraction of the progenitor disk. 
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Figure 16. Star-formation history and gas consumption for the G3G2 (top) and G3G1 (bottom) unequal mass interactions when merged 
on orbits that have varying amounts of angular momentum, as set by the pericentric distance R pcr i, which is listed in the top of each 
column in units of kpc. All merger orbits are inclined by 30° and nearly parabolic (eccentricity 0.95). The fiducial encounters all assumed 
Rperi = 13.6 kpc when the primary was G3. The burst efficiency is listed in the bottom of each column. All mergers use the nlmed 
feedback model. 



5.1 Variations in Orbit 

As outlined previously (Jl] and £|4.3[) . the enhancement of 
star formation during galaxy interactions originates from 
tidal forces that attend the merger. In particular, the close 
passage between the interacting galaxies produces bar-like 
structures in the stellar disk that torque gas into the galaxy 
center. Figure [8] shows that this occurs in both the primary 
and the satellite. The nearly prograde mergers that we have 
tracked up to this point produce a relatively strong bar ow- 
ing to resonances between the orbits of stars in the primary 
disk and the passage of the satellite galaxy. However, alter- 
native orbits and orientations of the interaction will affect 
the strength of the tidal force as well as the resonances ex- 
cited within the disk, and therefore the resulting starburst. 

To determine how the star-formation history depends 
upon the orbital angular momentum and disk orientation of 
the merging galaxies we have run two series of interactions 
using the G3G2 and G3G1 mergers. In the first series, we 
systematically change the orbital angular momentum which, 
in practice, is modulated by the pericentric distance R pcr i. 



All other parameters are identical to the fiducial mergers. 
The star-formation histories and gas consumption fractions 
resulting from this set of mergers are shown in Figure 1161 
For both the G3G2 and G3G1 interactions the maximum 
gas consumption occurs when R pcr i = 6.8 kpc, or roughly 
(6.8/2.85 w) 2.4 times the stellar scale radius of G3. The 
increasing gas consumption for more direct interactions sup- 
ports the connection between tidal forces and merger-driven 
star formation. We also note that the two interactions with 
Rpcri less than 6.8 kpc consume slightly less gas than orbits 
with larger R pcr i. This trend is likely a result of the disrup- 
tive nature of nearly head-on collisions as well the difficulty 
of capturing shock-induced star formation within SPH (al- 
though recent star formatio n model s have be en form ulated 
to address this shortcoming; iBarnesi (see, e.g.. |2004T0 . 

The second series of runs systematically alters the satel- 
lite orientation. While the fiducial set of mergers assumes 
that the orbital plane was inclined by 30° with respect to 
the disk of the primary galaxy, this series varied the inclina- 
tion. A range of mergers are performed in which the incli- 
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Figure 17. Star-formation history and gas consumption for the G3G2 (top) and G3G1 (bottom) unequal mass interactions when merged 
on orbits that have varying satellite orientations, from 0° (a perfectly prograde encounter) to 180° (a perfectly retrograde encounter). 
The other parameters of the merging orbits are similar to the fiducial case, i.e., they are nearly parabolic (eccentricity 0.95) and have 
moderate angular momentum (R pcr i = 13.6 kpc). Note that the fiducial set of mergers were inclined by 30°. The burst efficiency is listed 
in the bottom of each column. All mergers use the nlmed feedback model. 



nation is prograde (0°), polar (90°), and retrograde (180°), 
as well as several in between. Figure [17] shows the resulting 
star formation and gas consumption for these runs. As ex- 
pected, mergers in which the orbit of the secondary is aligned 
with the primary disk, i.e., prograde orientations, produce 
the strongest tidal responses in the disk and therefore the 
largest bursts of star formation. 



It is likely that similar results would follow from other 
orbital explorations, such as the initial separation between 
the two galaxies, as well as the energy of the initial orbit. In- 
terestingly, the largest starbu rst event in the series of major 
merger interaction studied bv lCoxl (I2004T ) was a zero net an- 
gular momentum orbit. In short, we expect any factor that 
efficiently strips gas of its angular momentum, including col- 
lisions, or increasing the tidal coupling between the satellite 
orbital motion and stars in the disk of the primary will lead 
to the most substantial bursts of star formation. 



5.2 Variations in the Extent of Gaseous Disk 

Moti vated by observations (see, e.g., lBroeils fc van Woerdenl 
1994), the galaxy models that form the basis for this work 
assume that the distribution of gas in the disk is more ex- 
tended than the stars. Specifically, the gas has an expo- 
nential scale radius a — 3 times that of the stellar disk. 
As a result, there is a significant quantity of gas at large 
radii. This extended gas does not significantly contribute to 
the nuclear starburst. Instead, the merger remnant contains 
a large quant i ty of gas in bot h a sh ock-heated hot phase 
jStruckl Il997l ; ICox et all |2004 120061 ) as well as in a cold 
star-forming disk. 

In order to discern how the assumed distribution of gas 
influences the burst efficiency, we have run a number of tests 
that alter the distribution of cold gas in the primary G3 
galaxy. Figure [18] shows an example of our typical result, 
namely, the extended gas distribution leads to an increased 
estimate of the burst efficiency. For the case shown in Fig- 
ure !18l the fiducial G3 galaxy model has been altered to have 
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Figure 18. Star formation and gas consumption for two G3G2 
mergers. The case on the left is similar to the fiducial interaction 
except the primary G3 galaxy model now has a = 1, where a is 
the multiplicative factor that sets the scale radius of the gaseous 
disk with respect to the stellar disk (see [[3J. The case on the 
right is also similar to the fiducial interaction except the primary 
G3 galaxy model has a = 3 and all gas beyond 20 kpc has been 
removed from the initial disk. The burst efficiency is indicated in 
the bottom panel for each case. Both mergers are run using the 
nOmed feedback model. 



a = 1, while keeping the remaining parameters, including 
the total gas mass, unchanged. In short, the vigorous star 
formation in the isolated primary consumes a larger fraction 
of the available gas and results in a lower burst efficiency. 

In additional experiments we merged models with the 
fiducial G3 distribution of gas (a = 3), however all gas 
beyond a cutoff radius (R cu toff) is removed. In essence, 
these galaxy models have much less gas, yet the inner, 
star-forming gas distribution is unchanged. As shown in 
Figure 1181 the resulting star-formation histories are un- 
changed. Even though the fractional gas consumption in- 
creases markedly (compared to the second column in Fig- 
ure [14]), owing to a corresponding increase in the gas con- 
sumption of the isolated system, the resulting burst efficien- 
cies are nearly identical. These experiments indicate that 
the burst efficiency is predominantly a function of the den- 
sity distribution of gas in the progenitor systems rather than 
purely the spatial distribution of gas. 



5.3 Variations in Bulge-to-Disk Ratio 

One of the most signific ant results from the work of 

iMihos fc HernquistJ |l994al . [ 1996) was the discovery that 
the internal structure of the primary galaxy can strongly 
influence the efficiency of the starburst. Specifically, these 
works showed that a massive stellar bulge stabilizes the disk 
against tidal perturbations and suppress strong inflows of 
gas that lead to starbursts. In the case of minor mergers, 
the presence of a bulge may eliminate the merger-driven 
starburst completely, while during a major merger the bulge 
may simply delay the starburst until the final coalescence of 
the two galaxies. 
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Figure 19. Similar to Figure [T5l except here the burst efficiency 
is shown for three series of runs with different values of the bulge— 
to-disk ratio (B/D) in the primary G3 galaxy. 



The absence of any discernible burst of star-formation 
for the G3G1 and G3G0 minor mergers presented in ^sup- 
port the notion that the presence of a stellar bulge can 
suppress merger-driven star formation when the mass ratio 
is large. Furthermore, this burst suppression persists even 
though the bulge adopted for the fiducial G3 model is only 
20% the mass of the stellar disk, i.e., the bulge-to-disk ratio 
(B/D) is 0.2, or about 50% smaller than in the models of 
Mihos & Hernquist where B/D— 0.33. In order to provide 
a more direct comparison to their work and further assess 
the effects of bulge mass on the merger-driven starburst, 
we have run a number of additional interactions where the 
bulge mass of the fiducial G3 model is altered such that 
there is no bulge at all (B/D= 0.0), or where the bulge mass 
is increased by a factor of 2.5 (B/D= 0.5). All other param- 
eters, including the fiducial merging orbits and the satellite 
galaxies (G2, Gl, and GO), are left unchanged. 

Figure [19] shows the burst efficiency as a function of 
merger mass ratio for our series of mergers in which the 
bulge mass of the G3 primary is varied. As expected, there is 
a systematic correlation for large mass ratio mergers to have 
smaller burst efficiencies when the bulge mass increases. For 
example, the bulgeless G3G1 (5.8:1) merger has a burst ef- 
ficiency of 0.16, which is three times larger than the burst 
efficiency for the model with the most massive bulge (B/D= 
0.5). We also note that the burst efficiency is insensitive to 
the bulge mass when the merger mass rati o is near unity, in 
agreement with I Mihos fc HernquistJ (|l996T ) . 

Using Eq. [5] the best fit to each series of B/D mergers 
is determined and overplotted in Figure All fits have the 
identical value of em, 0.49, a result of the constant burst 
efficiency during major mergers. The value of 7, which sets 
the mass ratio dependence, is 0.61, 0.74, and 1.02 for the se- 
ries with B/D= 0.0, 0.2 (the fiducial), and 0.5, respectively. 
It should be pointed out that the best fit for the fiducial se- 
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ries is slightly different than that reported in £|4.4. II because 
here only the G3Gx series is analyzed. 

Even though the qualitative relationship between 
bulge mass and burst effi ciency is in agreement with 
iMihos &: Hernquisd l|l994al ). their bulgeless minor merger 
produced a significant burst of star formation. The result 
is a burst efficiency (~ 0.7) that is much larger than our 
bulgeless series. Actually, this value is larger than any of 
our interactions, major or minor. This discrepancy was also 
noted in i]4.4.1l and is likely a result of three differences 
between their modeling and ours. 

First, the feedback model and the newer entropy- 
conserving version of SPH employed here both re sult in less 
inten se episodes of merger-driven star formation l|Cox et al.l 
2006). Second, the models employed by Mihos & Hernquist 
adopt a circular orbit for the satellite, which increases the 
tidal coupling to stars in the disk and produces a larger re- 
sponse. Lastly, the large burst efficiency found by Mihos & 
Hernquist is inflated by the very inefficient star formation 
assumed to occur in their quiescent disk. Such levels of star 
formation ap pear to be insuf ficient to match more recent 
observations l|Kennicuttlll998l ). 

While these three arguments outline why our burst ef- 
ficiencies are more modest — and we believe, more accurate 
— than prior calculations, the previous sections of this pa- 
per suggest that a number of the parameters for our fiducial 
encounters are sub-optimal at producing the largest merger- 
driven star formation event (and hence burst efficiency). We 
therefore performed a small number of additional simula- 
tions using a bulgeless version of our fiducial G3G1 merger, 
only we placed the satellite Gl on a co-planar (0°), close 
passage (R pcr i = 6.8 kpc) orbit in order to maximize the 
merger-induced starburst (this orbit was motivated by the 
results of i|5.1|l . 

The star-formation history during the interaction with 
the maximum burst efficiency, as determined from our small 
G3G1 parameter search, is shown in Figure 1201 We have 
also shown corresponding interactions when the modified 
bulgeless G3 model has bulge-to-disk ratios of 0.2, similar 
to the fiducial G3 model, and 0.5. Even with the large degree 
of parameter manipulation the maximum burst efficiency for 
the bulgeless run is only 0.19. While this is ~ 20% larger 
than the the bulgeless G3G1 in the fiducial encounter, and 
about ~ 4 times the fidu cial G3G1 merger, i t is sti ll far 
below the ~ 0.70 found bv lMihos fc Hernquistl i|l994aT ). 

5.4 Variations in Gas Fraction 

In our final set of additional merger simulations, we vary the 
gas fraction of our fiducial G3 galaxy model. To this end, the 
total mass in the disk is kept fixed, yet the distribution of 
mass, i.e. the amount in the gaseous versus the stellar disk, 
is varied for each model. All other parameters, including 
the interaction orbit and the satellite galaxies (G2, Gl, and 
GO) remain unchanged from their fiducial values. Figure [21] 
presents the star formation history, the gas consumption, 
and the burst efficiency for the fiducial G3G2 interaction 
plus two additional mergers when the G3 primary has a 
larger gas fraction, /, defined as the mass of the gas disk 
divided by the total (gas plus stars) disk mass. 

Figure [21] exemplifies the generic outcome of all inter- 
actions that test the gas fraction of the primary, namely, 
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Figure 20. Star formation during unequal mass interactions as 
a function of primary bulge mass. The interaction is similar to 
the fiducial G3G1, except the orbit is co-planar (0°) and has less 
angular momentum (see text) in order to produce the maximum 
starburst. All runs are performed with the the nOmed feedback 
model. 



increasing the primary's gas fraction decreases the burst effi- 
ciency. In effect, the result is similar to that presented in t|5.2l 
for the distribution of gas and follows from a similar cause 
— the increased gas consumption in the isolated primary. In 
both scenarios, the initial disk converts a large fraction its 
gas into stars regardless of whether the interaction occurs 
or not. 



6 DISCUSSION 

In this paper we perform a series of numerical simulations 
that follow the interaction and merger of binary galaxies 
with various mass ratios. Our analysis quantifies the star- 
bursts that result from the tidal forces that attend the 
merger. As expected, mergers between galaxies with nearly 
equal mass generate the largest tidal forces, and therefore 
produce the most intense bursts of star formation. Mergers 
between galaxies with a large mass ratio produce relatively 
little, if any, enhancement in star formation above quiescent 
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Figure 21. Star formation history and gas consumption for the 
fiducial G3G2 merger and two similar interactions where the pri- 
mary G3 galaxy has a larger gas fraction. As with Figure the 
difference between the merger and isolated evolution is shaded for 
clarity. The burst efficiency is displayed in the lower-right of the 
bottom panel. We note that the simulation here were run with 
the n2med feedback model while the mergers shown in Figure IT4l 
employed the nOmed model. 



evolution, yielding a correlation between merger-driven star 
formation and decreasing merging galaxy mass ratio. 

To quantify the relationship between merging galaxy 
mass ratio and merger-driven starburst, we introduce the 
burst efficiency (Eq. [5} as the fraction of gas that is con- 
verted into stars during the interaction which does not occur 
if the galaxies are evolved in isolation. The burst efficiency is 
demonstrated to be robust to uncertainties in the feedback 
parameterization unlike the strength or duration of the star- 
burst. 

While the burst efficiency resulting from collisions be- 
tween galaxies of equivalent mass is relatively insensitive 
to the details of the merging event such as the orbit, the 
galax y orientation, and the properties of the merging g alax- 
ies (|Mihos fc Hernquisull996l ; ISpringelioOol ; ICoxlbool ) . this 
is not the case when the participating galaxies are unequal 
in mass. 

By performing a large number of additional merger sim- 
ulations, we have quantified the effects of merging orbit and 
orientation, as well as properties of the progenitor disk. 
We find that close passage, co-planar orbits produce the 
most significant bursts of star formation, consistent with 
the expectations of the tidally induced origin for the star- 
bursts. We also find that the structure of the progenitor disk 
strongly influences the merger-driven star formation. In par- 
ticular, the presence of a centrally-concentrated stellar bulge 
stabilizes the disk and suppresses merger-driven star forma- 
tion. The distribution and mass of the gaseous disk also 
influences the starburst. In general, increasing the amount 
of gas at densities above p cr it, the threshold density for star 
formation to commence, decreases the burst efficiency. 

These results lead us to conclude the following two 
facts about merger-driven starburst during an unequal mass 



galaxy merger: (1) significant starbursts occur for only spe- 
cialized scenarios, e.g., close passage, co-planar orientations, 
when the primary disk does not contain a bulge, and (2) 
even in this extreme scenario, the burst efficiency is still 
only < 0.25, i.e., a single unequal mass merger does not 
convert a large fraction of gas into stars. 



6.1 Comparison to Previous Simulations 

Our work is closely related to, and consistent with, a number 
of prior studies of the gas dynamics durin g the interaction 
and merger of unequal mass galaxies (e.g., Hernquistlll989l ; 
iMihos fc HernquistJ|l994al ; Irlernquist fc Mihodll995f rWhile 
all work performed to date demonstrates that minor mergers 
can induce radial inflows of gas that result in periods of en- 
hanced star formation, our models typically result in smaller 
burst efficiencies (see H4.3.1I and ij4.4.1|) owing to three key 

differences in the studies. 

First, as shown in lCox et al. (2006) using simulations of 
major mergers, the more complex treatment of the interstel- 
lar medium produces more hot gas, and suppresses star for- 
mation. In particular, relaxing the isothermal gas assump- 
tion commonly employed in the work of Mihos fc Hernquist, 
inclu ding the "conservative-entropy" (|Springel fc Hernquistj 
I2002T) version of S PH, and the more efficient feedback models 
of I Cox et al.l l|2006h all serve to suppress the merger-induced 
starburst, and lower the burst efficiency. 

Second, the se t of simulatio ns employed by 
IMihos fc Hernquistj (fl994a) and Irlernquist fc Mihod 



( 1995T ) followed satellite galaxies that were initially placed 
on a circular orbit. In contrast, we follow parabolic orbits 
which are motivated by cosmological expectations (see Q. 
These more energetic orbits lead to less direct coupling 
between the orbital angular momentum of the satellite and 
the disk of the primary, and therefore are less conducive to 
an intense inflow of disk gas. 

Lastly, as noted in £15.31 and in ICox et all (|2006t ). the 
star formation model employed by Mihos & Hern quist is less 
efficie nt than current observations suggest (e.g., iKennicvrtq 
1998). The primary result of this assumption is that the 
quiescent galaxy consumes much less gas when evolved in 
isolation and the merger burst efficiency is overestimated. 

While there are additional differences between the 
galaxy models employed by these previous studies, namely 
they use less massive and less concentrated dark matter ha- 
los, and the baryonic components are not as faithful a rep- 
resentation of observed galaxies in the local Universe, the 
tests performed in iJS] indicate that these differences play a 
secondary role to the three items outlined above. Even when 
we modified every possible parameter to maximize the burst 
efficiency, our value was still only one-third as large as the 
previous results. 



6.2 Comparison to Observations 

While observational studies have already established 
a clear link between star formation and galaxy inter- 
actions (as measured by clo s e pairs or morp h ology , 
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has it been possible to specifically address whether or not 
this correlation holds when the close galax y pairs have 
rather disparate luminosities. In particular, IWoods et al.l 
(2006) reported no correlation between star formation 
(as measured by Ha) and close galaxy pairs in the CfA2 
Redshift Survey when the magnitude difference is greater 
than 2, and a more recent study using a larger sample from 
SDSS found that the satellite galaxy indee d shows enhanced 
star formation while the primary did not l|Woods fc Gellerl 
120071) . 

Assuming that luminosity traces mass (and modulo any 
systematic affects owing to merger-driven star formation), 
the studies by Woods et al. imply that there is no obser- 
vational evidence for induced star formation in the primary 
galaxy when the merger mass ratio is greater than ~ 6:1, 
while the satellite is more likely to experience a starburst. 
Taken at face value, our results naturally recover this ob- 
served trend. Our merger models produce little, if any, glob- 
ally enhanced star formation for merger mass ratios below 
5:1, and when they do, it requires very specific circumstances 
(bulgeless primary, co-planar, close-passage orbit). It is also 
intriguing that Figure[8]hints that the satellite is much more 
susceptible to enhanced star-formation during the interac- 
tion, which also seems consistent with the observations. 

In contrast to the statistical studies, there are a mul- 
titude of observations that suggest that individual systems 
are currently undergoing minor merger induced episodes of 
star formation (see the list in Jill. To determine whether 
these specific galaxies are consistent with the statisti- 
cal studies requires more extensive modeling of individ- 
ual systems, a s has been performed in a few cases al- 
ready (see, e.gjMihos fc Bothunlll997l ; lLaine fc HeUeill 19991 : 



ISalo fc Laurikainenl 12000). It is possible that many of these 
systems have unique satellite orbits, or that the effects of 
multiple minor mergers which occur simultaneously is more 
dramatic than the binary mergers that we have followed 
here. Models of individual galaxies, and their direct compar- 
ison to observations across many wavebands may also yield 
important constraints to the star formation and feedback 
models. 

6.3 Implications for Dwarf Galaxies 

While the focus of our analysis has been on the global 
properties of merger-driven star formation, Figure [5] clearly 
demonstrates that the star formation history of the satellite 
galaxy can be enhanced far more dramatically than that of 
the primary galaxy (which dominates the global star forma- 
tion). Such a scenario appears to be observed in th e large 
sample of SDSS galaxies studied by I Woods fc Geller! (|2007m 
and also in the nearby Universe, e.g., both M82 and NGC 
3077 are currently experiencing periods of intense star for- 
mation after recen t close passages t o M81 (|Yun et al.lll994l : 
I Walter et al]|2002j : lott et al.l 120031 ) . These tidally-induced 
episodes of star formation may also have implications for 
the detection of satellite g alaxies and the infe rred cosmolog- 
ical merger rate (see, e.g.. iBerrier et aI1l2006l ). 

6.4 The Mass— Dependence of Star Formation 

One interesting feature present in Figure [15] is the system- 
atic dependence of ei : i, the burst efficiency for an equal- 



mass major merger, on primary galaxy mass. Specifically, 
the burst efficiency is 0.46 for the G3G3 major merger and 
steadily increases to 0.61 for the G0G0 major merger. 

The increasing burst efficiency era with decreasing sys- 
tem mass is a direct byproduct of systematic changes in the 
merger-induced star formation compared to that in the iso- 
lated disks. One possibility for this trend is the systematic 
increase in gas fraction with decreasing galaxy mass that is 
assumed for our galaxy models. However, the results of >15.4I 
suggest that increasing the gas fraction actually decreases 
the burst efficiency. 

A more likely scenario is that the density-dependent 
description of star formation, including the explicit density 
threshold p ZT i t for star formation to commence, is producing 
the large variation in burst efficiency with mass. In par- 
ticular, Figure [2] shows that the quiescent star formation 
spans three orders of magnitude from G3 (~ 1 Mq yr _1 ) to 
GO (~ 1O _3 M0 yr _1 ), while the peak star formation dur- 
ing the major mergers are much more comparable (G3G3 
~ 25M Q yr _1 and G0G0 ~ 3M Q yr" 1 ). The significant in- 
crease in merger-induced star formation compared to qui- 
escent levels was noted in Q4.3\ which determined that the 
G3G3 merger enhances star formation by a mere order of 
magnitude while the G0G0 merger enhances star formation 
by a whopping three orders of magnitude. 

More work is required to determine the precise nature of 
the mass-dependence of star formation and how this trend 
depends upon uncertainties associated with our implemen- 
tation of star formation. Moreover, a systematic compari- 
son to the observed relationship between specific star forma- 
tion rate and stellar mass, i.e. downsizing l|Cowie et~al ]|l996l : 
iBrinchmann fc Ellisll200ol ; iKauffmann et alj|2003l ), may elu- 
cidate the physical mechanisms responsible. 

6.5 Input for Future Studies 

Quantifying the merger-driven star formation as a func- 
tion of merger mass ratio is useful for a variety of further 
studies. In particular, semi-analytic models of galaxy for- 
mation (SAMs) often find that merger-driven star forma- 
tion is necessary to reproduce the lumino sity function and 
numb er counts of Lyman-break ga laxies jS omcrvillc et ail 
200lT) and sub-mil limeter galaxies (|Guiderdoni et aT, 19981 : 
Baug het al.l 120051). To this end, we have followed the 
methodology of ISomerville et"al] (|200ll ) and introduced the 
burst efficiency (see t|4.4.1l and Eq. [5} in order to param- 
eterize the star formation enhancement that occurs during 
galaxy interactions as a function of participant mass ratio. 

In Table [3] we list the best fit parameters to Eq[5] for 
all models explored in this paper, including the fiducial se- 
ries of runs ( H4.4.ip . and all additional models (orbits in 
E{5TT] bulge-to-disk ratios in i]5.3l and gas fractions of the 
progenitor disk in £|5.4[1 that use the G3 primary galaxy. 
A quick inspection of Table [3] indicates that the parameter 
ei : i, which normalizes Eq.[5]to the burst efficiency of a ma- 
jor merger, changes relatively little across this wide range 
of tests, although it may depend on galaxy mass (see the 
discussion of i]6.4[) . On the other hand, the parameter 7, 
which determines the mass ratio scaling, has a significant 
dependence on orbital inclination and bulge-to-disk ratio, 
a mild dependence on the orbital angular momentum, and 
a negligible dependence upon gas fraction. 
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Table 3. Compilation of best fit burst efficiency (see £|4.4,ll and 
Eq. parameters for various sets of simulations employed in this 
paper. All simulations analyzed for these fits used the n2med 
feedback model. 
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0.49 
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EH 
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EH 


0.50 


0.74 


B/D= 


= 0.2 


E3 


0.50 


1.02 


B/D= 


= 0.5 


15^1 


0.50 


0.74 


/= 0.20 


EH 


0.49 


0.74 


/= 0.50 


EH 


0.44 


0.72 


/= 0.78 



While the burst efficiency is a useful quantification of 
the star-formation induced by a galaxy merger, a better 
understanding of the star-formation timescale is necessary 
to completely describe merger-induced star formation. Un- 
fortunately, as shown in ^4.3.31 uncertainties in the feedback 
model do not permit an unambiguous characterization of the 
star-formation timescale. There is hope, however, that bet- 
ter constraints can be placed on the feedback model through 
more detailed modeling of individual systems systems or by 
comparison to the observed distribution of star formation 
rates. 

Finally, we emphasize that the accretion events followed 
in this study drive structural evolution which is also an in- 
teresting and relevant input parameter to future studies of 
galaxy formation and evolution. In particular, the morphol- 
ogy of the merger remnants (see Figs. [3H6} suggests that 
the stellar and gaseous components react differently to the 
dynamical perturbation and that the remnant appears to 
be a systematically earlier Hubble type than the original 
disk. This is also consistent with the surface density profile 
shown in Figure [5] which shows an excess of mass at small 
radii indicative of the formation of a stellar bulge. 



6.6 Other Considerations 

As a final comment, we note that the present simulations 
lack several physical processes that may play a role in the 
results we have presented. First, we have not included the 
recycling of gas from stellar winds and supernovae, an omis- 
sion that would increase the gas fraction as a function of 
time. Another process that may increase the gas content on 
long timescales is the accretion of gas from the cosmological 
growth of structure. While we eventually intend to incor- 



porate these effects into our modeli ng (as has been done 
by several authors already, see, e.g, Tornatore et al.l |2004| ; 
IScannapieco et ail 120051 ; IStinson et al.ll2006l ). the relatively 
small difference (~ 15%) in the burst efficiency for runs of 
different initial gas fractions (see i)5.4|l supports the notion 
that this omission has a relatively minor effect on the star 
formation that results from a single merger. However, for 
the long-term evolution, and for the absolute star-formation 
rate, these effects are likely to be important. 

It should also be pointed out that our galaxy models are 
calibrated to match low redshift observations of disk galax- 
ies. While we have surveyed a small portion of the parameter 
space that might correspond to properties of disks at higher 
redshift, e.g., higher gas fractions and more compact ini- 
tial disks, and therefore have some indication of the changes 
that may occur, future studies will need to investigate these 
dependencies in more depth. 

We have also not included accreting black holes and 
their associated feedback into the models presented here. 
Recent work has shown that these processes may play a sig- 
nificant role dur ing major mergers, r esulting in the forma- 
tion of quasars ([Hopkins et al.l 120060 . and leading to rem- 
nants that reside on the Mbh — o relation (|Pi Matteo et al.l 
120051 ; ISpringel et al.l l2005bl; IRobertson et all I2006T ) and 
have the colors (ISpringel et al.l l2~005al) an d other proper- 
ties l|Cox et al . 2006; IRobertson et alJuOOfj ) appropriate for 
present day elliptical galaxies. While the minor interactions 
we discuss here are a likely fueling mechanism for many 
form s of nuclear activity (see, e.g., iHopkins fc Hernquistl 
|2006j), because significant black hole growth requires large 
quantities of gas (approximately the entire content of the ini- 
tial disks) to be driven to the galactic center, it is likely that 
black hole growth and feedback play a relatively minor part 
during most of the unequal mass interactions we follow here. 
Of course, the existence of a correlation b etween the black 
hole mass and the bulge mass (see, e.g., iMagorrian et al] 
1998) suggests that at the very minimum the black hole im- 
parts enough feedback to regulate its own growth. 
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